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ABSTRACT 

Neutrino radiation takes a major role in the momentum, heat, and lepton transports in proto- 
neutron stars (PNSs). These diffusive processes affect the growth of magnetorotational instability 
(MRI) in PNSs. We perform a local linear analysis for the axisymmctric and nonaxisymmetric 
MPJ including the effects of neutrino transports and ohmic dissipation. We find that the MRI 
can grow even in the multi-diffusive situations that are realized in neutrino loaded PNSs. When 
the toroidal magnetic component dominates over the poloidal one, nonaxisymmetric MRI modes 
grow much faster than axisymmetric modes. These results suggest the importance of the nonax- 
isymmetric MRI in PNSs. Thus the understandings of three-dimensional nonlinear evolutions of 
the MRI are necessary to reveal the explosion mechanism of core-collapse supernovae. 

Subject headings: instabilities — MHD — stars : magnetic fields — stars : neutron 

1. INTRODUCTION 

Strong magnetic fields can play an important role on the dynamics of core-collapse supernovae. Akiyama 
et al. (2003) pointed out the possibility of the field amplification by the magnetorotational instability (MRI) 
in supernova core. The shapes of shock wave and the neutrinosphere are modified aspherically by the 
coupling effects of the rotation and magnetic fields, and anisotropic neutrino radiations could lead to jet- 
like explosions (Kotake et al.2003, 2004). Magnetic pressure driven explosions are also investigated in the 
context of magneto-rotational core collapse (Yamada & Sawai 2004; Mizuno et al. 2004; Takiwaki et al. 
2005; Ardeljan et al. 2005; Sawai et al. 2005; Moissenko et al. 2006). In these studies, strong magnetic fields 
(> 10 15 G) are generated by the field compression, field wrapping, and the MRI, which lead to the prompt 
explosions of core collapse supernovae. Alternative effects of magnetic fields are proposed by Thompson et 
al.(2005). They suggest that the turbulent viscosity sustained by the MRI converts the rotational energies 
of supernova core into the thermal energies, and thus the MRI will affect the dynamics of core collapse 
supernovae. However, these studies mainly focus on the dynamical effects of magnetic fields at the outside 
of PNSs. Therefore, magnetic phenomena inside PNSs and its influence on the core collapse dynamics have 
not been understood yet. 

There are a few studies of magnetic effects at the interiors of PNSs. Miralles et al. (2002) show 
that the convective fluid motions in PNSs are prevented by strong magnetic fieds. Numerical simulations 
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suggest that the magnetic bubbles rise up by buoyancy and carry neutrino-rich materials to the neutron star 
surface (Wilson et al. 2005). This could increase the neutrino luminosity sufficictly to achieve a successful 
neutrino driven explosion. Socrates et al.(2005) show that the neutrino bubble instability grows at the 
inside of magnetized PNSs. Nonlinear evolution of the neutrino bubble instability will generate large density 
fluctuations. The radiation preferentially leaks out of the rarefied regions, which enhances the outgoing flux. 
They also propose that these global asymmetry in the neutino emission would cause the pulsar kicks. 

Masada et al. (2006; hearafter Paper I) investigate the growth of the nonaxisymmetric MM and show 
that it can grow even in the stably stratified envelope of PNSs. Our results suggest that the nonlinear 
evolution of MRI amplifies the magnetic fields and drives MHD turbulence in PNSs. The layers below the 
ncutrinosphcre are thought to be convectively stable. If buoyant motions are excited by the MRI, it could 
lead to enhancement of the neutrino luminosity. However, in Paper I, the effects of neutrino radiations 
on the growth of MRI are neglected. Below the neutrinosphere, electron-type neutrinos exchange energy 
with the matter via pair capture process and exchange momentum via elastic neutirno-nucleon scattering. 
These processes dominate the other transport processes by electrons and/or photons in PNSs, and thus the 
coefficients of the momentum, heat, and lepton transports are determined by neutrino radiative processes. 
The importance of these transports caused by the neutrino radiation have been suggested already in the 
previous studies of the mixing processes in PNSs (Bruenn & Dinneva 1996; Mezzacappa et al. 1998; Miralles 
et al. 2000; Bruenn et al. 2004; Socrates et al. 2005). Therefore, we must study such neutrino effects to 
reveal the growth of MRI in more realistic, neutrino loaded PNSs. 

In this paper, we analyze the diffusive effects caused by neutrino radiations on the growth of MRI and 
discuss the natures of MRI in PNSs. We consider four kinds of diffusive processes; the momentum, heat 
and chemical diffusions induced by the neutrino radiation and the ohmic dissipation of hot nuclear matters. 
As the mean free path of neutrino is short compared to the density and temperature length scales below 
the neutrinosphere, we can use diffusion approximation for the neutrino transports (Socrates et al. 2005). 
Another important feature of our analysis is to consider the vertical magnetic field as well as the toroidal 
one. For this case, both axisymmetric and nonaxisymmetric modes can be unstable for the MRI. In previous 
works, axisymmetric and nonaxisymmetric MRI are investigated separately (Acheson 1978; Balbus & Hawley 
1994; Menou et al. 2004). We consider both axisymmetric and nonaxisymmetric stability and discuss their 
difference in the growth rate and stability conditions. 

The paper is organized as follows. In § 2 we obtain the dispersion equation including diffusive processes, 
which determines the stability of magnetized PNSs. Necessary conditions for the stability are derived by an 
analytical approach. In § 3 we discuss general properties of the dispersion relation focusing on the effects 
of quadruple diffusivities. In § 4 we apply our results to the interior of magnetized PNSs, and examine the 
effects of neutrino radiations on the growth of MRI. We also examine the difference between axisymmetric 
and nonaxisymmetric modes of MRI. The possibility of the other instabilities is discussed in § 5. Finally, we 
summarize our main findings in § 6. 

2. LOCAL DISPERSION EQUATION 

2.1. Physical assumptions and basic equations 

We examine the stability of differentially rotating, magnetized PNSs using a linear perturbation theory. 
In this paper, we focus on the effects of the neutrino radiation for MHD instabilities. The weak interaction 
between nuclear matters and neutrinos is described by using various opacities. We simply assume that 
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the energy and momentum are exchanged between electron-type neutrinos and nuclear matters via pair 
capture process and elastic neutrino-nucleon scattering. In addition, neutrino opacities are assumed to be 
proportional to the square of the neutrino energy E, 

k v = k (E/E ) 2 , (1) 

where Ko and Eq are constants and independent of composition (e.g., Socrates et al.2005). In this paper, 
opacities arc defined as the inverse of mean free paths. 

Neutrinos are optically thick and both the thermal and chemical equilibrium would be held below the 
neutrinosphere. Thus, the governing equations are the following, 

^+V-(pu) = 0, (2) 
du . , _ 1„ 1 _ 2 \ . 1 



m + (u • V> = --V (P + + -(B.y )B + g + v^u , (3) 

— = V x (u x B) + r]V 2 B , (4) 
n(^+u-VY L ) =-V-F L , (5) 

nT(^+u- Vs) + n M „ e + u • VF L ) = -V • F . (6) 

These are the continuity equation, equation of motion, induction equation, conservations of lepton fraction 
and energy, respectively. Here p is the fluid density, u is the flow velocity, P is the total pressure of nuclear 
matters and radiations, B is the magnetic field, g is the gravitational field, n is the baryon number density, 
s is the entropy per baryon, Yj, is the lepton fraction, p Ve is the neutrino chemical potential, and F and Fl 
are the radiative energy flux and lepton flux, respectively. Neutrino viscosity and magnetic diffusivity are 
represented by v and rj. The effect of the bulk viscosity is neglected here. 

We employ the equilibrium diffusion approximation, and thus our analysis is restricted to the region 
below the neutrinosphere. Using the energy flux of the electron neutrino and anti-neutrino, F Vt and F^ , 
the radiative energy flux F is given by 

F = F Ve + F v - e = -a T VT - a L VY L , (7) 

where 



go k l T = go ^ c (dn v \ 

1 6 Ko h 3 c 2 3 ' L 6 Ko ft 3 c 2 tt 2 \dY L )p,T y > 



(Bludman & Van Riper 1978). The lepton flux is given by 

F L = -b T VT - b L VY L , (9) 

where 

^-'^-■sabMml,- <10) 

Using the simple parameterization of opacity defined by equation (1), the temperature gradient cannot 

contribute to the lepton flux, that is br = (Socrates et al. 2005). Using equations (7) and (9), we can 
rewrite the conservations of lepton fraction and energy; 

n(^+u-VY L ) =b L V 2 Y L , (11) 
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nT^+u- Vs) = a T V 2 T . (12) 

We ignore the spatial dependence ol the diffusion coefficients ar, o-l-, and 6l. This is appropriate for a local 
linear analysis. The resistive and viscous dissipation terms are also neglected in the energy equation, because 
they are higher order terms. 



2.2. Linearized Equations 

Following the analysis performed in Paper I, we consider the Eulerian perturbations (denoted by a 
prefix S) with the WKB spatial and temporal dependence, <5 oc exp{i(fc r r + ni(f) + k z z — at)}. Local WKB 
analysis is irrelevant in the cases that the growth rate of nonaxisymmctric unstable modes is smaller than 
the differential rotation rate. However, we now focus mainly on the MRI in the cases with strong velocity 
shear for which the growth rate is comparable to the angular velocity and adopt the simple WKB analysis 
in this paper. We use the cylindrical coordinates (r, <j> 1 z) and consider a PNS rotating with the angular 
velocity f2(r), and its magnetic fields being B = (0, B^, B z ). We assume magnetic fields are uniform locally 
for simplicity. At the unperturbed state, the PNS is in the magnetohydrostatic equilibrium, that is, 

B l IdP , IdP 

Z^Z = -9r + m\ zl r- = -g x , (13) 



Airpr p dr p dz 

in the radial and vertical directions, respectively. Written out in component form and the largest terms 
retained, equations (2)-(4) become to linear order: 



k r 5u r + k z 8u z = , (14) 

(15) 



i{k-B)„ ik r ik r (B ■ Sb) B^Sbs So. 

iw v Su r + 2CISuj> + -Sb r = —SP H ^ '- + J g - — (g r - rf2 2 ) 

47T/3 p 47T/9 2npr p 



lUvSuj, 8u r + — Sbj, + -f- = , (16) 

r or Anp Aitpr 

. . i{k-B) £u ik Z£n ik z (B-Sb) 5p 

iloJuz + —. Sb z = — SP + \ '- - —g z , 17 

4irp p 4irp p 

u) v Sb r + (k ■ B)Su r = , (18) 

d^l B 

iuj„Sb,f, + i(k ■ B^Sua, + r— — Sb r H -Su r = , (19) 

or r 

Lu n Sb z + (k ■ B)5u z = , (20) 



where 



u = a — mfl , lo u = oj + ivk 2 , uj v = lj + ir\k 2 , 

k 2 = k 2 r +k 2 , k = (k r , m/r, k z ) , Sb = (Sb r , Sb^, 5b z ) . 

Here we use the local approximation (k 3> l/r,m/r). We also adopt the Boussinesq approximation because 
the typical speeds considered in our analysis are much smaller than the sound speed. The relation SP + 
B^Sb^/An — is used in equation (16) by virture of the local approximation (Acheson 1978). 

To close linearized equations, we require another equation obtained from the thermodynamic processes. 
The chemical equilibrium is realized in PNSs, and thus the density is a function of pressure, temperature, 
and lepton fraction. In the Boussinesq approximation, pressure perturbations arc negligible because fluid 
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elements are assumed to be in the pressure equilibrium with their surroundings. Thus the density and 
entropy perturbations can be expressed in terms of the perturbations of temperature and lepton fraction: 



Sp 



p(a^ + (35Y L ) , 



ST 



5s = m B c p — + (5Y L , 



(21) 
(22) 



where mj is the baryon mass, a — — (dlnp/dlnT)p t Y L , = — (d In p / 8Yl) p,t are the coefficients of thermal 
and chemical expansion, and c p — (T/m B )(ds/dT)p t Y L is the specific heat at constant pressure. The sign 
of £ = {ds/dYL)p,T determines whether the leptonic gradient is stabilizing or destabilizing. 



Equations (11) and (12) can be linearized and written in the form 

-iu L 5Y L + 5u-VY L = 0, 
CT 

-iu T 6T - Su ■ V$ T = —b L k 2 SY L 
pc p 



where 



lot = oj + i\k , lol = lo + i^k , 

T fdT\ 

V $ T = (Vs-(VY L )= — VP-VT. 

m B c p \oP/s,y l 



(23) 
(24) 

(25) 
(26) 



We introduce the heat diffusivity \ = ax/pep and the chemical diffusivity £ = b^jn. Equation (26) represents 
the subadiabatic temperature gradient. Substituting equations (23) and (24) into equation (21), we obtain 
the linearized equation, 



iu^ujt = a^(Su ■ V$ T ) - (3[w + z(l + ^)xk 2 ](Su ■ VY L ) , 
p T 



where 



_ 3 p Ve fdlnp, 



/ 64n p Ue \ / d s/k B \ 
V 91nT )s,p\ dY L )t,p 



(27) 
(28) 



7T 2 k B T 

Assuming the nondegenerate neutrino and anti-neutrinos in degenerate nuclear matters, ip is considered to 
be sufficiently smaller than unity (Burrows & Lattimer 1986). Thus, equation (27) is rewritten as 



Sp ljl 

iuJTUJL — = a—(Su ■ V$t) — (3lot(Su ■ VYl) . 



(29) 



Eliminating perturbed quantities in eight linearized equations [eqs. (14) - (20) and (29)], we obtain the 
following dispersion equation: 



It y It y 



iV Te + i\ Le 

LOT 



2 2 2 




2lo 2 a 




LO D - K UJ R 


K 2 





-4fi 2 (fe • v A f 



2 V O 



)\k-v A ) 2\ Q 



k ■ v A 



(30) 



where 



J D 



v A = 



J v LO n - (k ■ v A ) 2 , 



\ = <4 - ( fe • v a? > 

^ = (0,V M ,V Az ), K 2 = ^(rW) 



B 
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A r 2 4lTpr 2 ' Az r2 Alrpr 2 > 



D = 



k 2 P T y ' J ' Le k 2 p 
' d k r d \ _ dlnQ 

^dr k z dz) ^ dlnr 
Here, va is the Alfven velocity, n is the epicyclic frequency, u>a is the azimuthal Alfven frequency, ujaz is the 
vertical Alfven frequency, Nt & is the effective thermal buoyancy frequency, and Nl s is the leptonic buoyancy 
frequency The shear parameter q denotes the rotational configuration of the systems. 

Since the dispersion equation (30) has a very complex form, it is difficult to treat for what it is. Therefore, 
we simplify it using reasonable approximations. When we focus on the MRI modes with u <~ k ■ va and 
adopt the weak magnetic field approximation £1 3> uja, equation (30) can be described as 



k 2 k 2 

T2 U D - 



Z2LN 2 i ZH N 2 



ui z D - k z lu z r - 4(fc • v A ) ^ = . (31) 



It is stressed that both axisymmetric and nonaxisymmetric MRI modes are included in this dispersion 
equation. In the diffusionlcss {v = il = x = £ = Q) and pure toroidal (B z = 0) limits, equation (31) is 
identical to the result of Paper I. Taking chemically homogenious (VYl = 0) and axisymmetric (to = 0) 
limits, it reduces to the dispersion equation derived by Menou et al.2004. 



2.3. Stability Criterion 

The nature of instabilities in PNSs has been considered by a number of authors (Brucnn & Dinneva 1996; 
Mezzacappa et al.1998; Miralles et al. 2000, 2002, 2004). We make a detailed comparison of our dispersion 
equation with previous studies and derive the stability criterion for the MRI in the neutrino loaded PNSs. 

Taking the diffusionlcss and hydrodynamic (B = 0) limits, equation (31) is identical to the result of 
Miralles et al.2004, in which the stability conditions for the differentialy rotating PNSs are obtained. In 
the non-rotating limit (il = 0), equation (31) reduces to the results of Miralles et al.2002, in which they 
derive general stability criteria for the convective instability, taking dissipativc processes into account such 
as neutrino transport, viscosity, and resistivity. Thus, the both semiconvective and neutron-finger unstable 
modes are obviously involved in our dispersion equation (Bruenn & Dinneva 1996; Miralles et al.2000). It 
is, however, noticed that we neglect the influence of the chemical inhomogencity on the heat diffusion and 
the thermal inhomogeneity on the chemical diffusion for simplicity. 

Although the linear and nonlinear growth of semiconvection and neutron-finger instability would be a 
quite important issue in the evolution process of PNSs, it is beyond the scope of this paper. We now have 
an interest with the effects of multiple diffuions on the growth of MRI. Hence, we show the general stability 
criterion for the MRI in multi-diffusive systems, which is derived from equation (31) with the same procedure 
as in Urpin 2006. When we focus on relatively larger wavelengths (A > A c e 2n y/vrj/VAz), the stability 
criterion for the MRI can be written as 

!L N^ + M + ^->0, (32) 
X £ dlnr 
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where the thermal buoyancy frequency Nt and the leptonic buoyancy frequency Nl arc defined by 

N*=-^-(VP)-Vy Tl Nl = ^(VP)-VY L . (33) 
pT p 

Equation (32) have a form similar to that derived in Urpin 2006 but with the leptonic gradient and the 
chemical diffusion. This criterion differs essentially from the standard stability criterion of the MR1 (Balbus 
& Hawley f 991, 1994). Considering the diffusive coefficients which are typical in PNSs (x,£,3> V- sec § 4.1), 
the stabilizing effect of the stratification is much reduced in our criterion. Here we notice that the viscous 
effect is not appeared in above criterion. However, the growth of MM can be suppressed by the viscous 
damping (see § 3.2). Thus we should consider not only the stability criterion but also the unstable growth 
rate when we discuss the characteristics of the MRI in multi-diffusive situations. 



3. GENERAL FEATURES OF THE DISPERSION EQUATION 

We solve the dispersion equation (31) numerically and examine the characteristics of MRI in multi- 
diffusive situations. Note that, within the local approximation, the dispersion equation (31) is valid in the 
case f2 > k v^ ^> &a- As described above, both axisymmetic and nonaxisymmetic MRI modes are included 
in our dispersion equation. General features of MRI in such complicated situations are not fully investigated 
in previous studies. Our analysis could be applicable to not only the PNSs but also the other astrophysical 
objects, such as stellar interiors and accretion disks. 

The effective thermal and leptonic buoyancy frequencies, NT e and Nlc denote the size of the restoring 
force due to the stable stratification. Qualitative effect of the leptonic buoyancy is the same as thermal 
one. Thus we consider only the thermal buoyancy frequency as the measure of stabilizing effects due to the 
stratification for a while. In PNSs, the toroidal magnetic fields would dominate over the poloidal one. In 
the following, we fix the ratio of the toroidal to poloidal Alfven frequencies as loaI^Az = 10 4 - The ratio of 
the angular velocity to the toroidal Alfven frequency is assumed to be CI/u)a = 100. These are about their 
typical values in rotating PNSs (see § 4). Since we have an interest to the systems with differential rotation, 
the shear parameter q = 1.5 is adopted in this section. Furthermore, we restrict our analysis to the modes 
with zero radial wavenumber (fc r = 0) which correspond to the fastest growing branch of the axisymmetric 
and nonaxisymmetric MRI. 



3.1. Diffusionless Limit 

We show the features of the dispersion relation in the diffusionless limit (^ = ?7 = x = <!; = 0). When 
the imaginary part of u> is positive, the amplitude of the perturbation will grow exponentially Figure 1 
shows a three-dimensional plot of the unstable growth rate as a function of the vertical and azimuthal 
wavenumbers. Stabilizing effects due to the stratification are ignored here (Nr e = 0). The growth rate is 
normalized by angular velocity J7. Both axisymmetric and nonaxisymmetric modes can be unstable for the 
MRI in this situation. For the case with q = 1.5, the maximum growth rate of axisymmetric modes is 0.7517 
at k z = k max = 0.97CI/Vaz- The cutoff wavenumber of axisymmetric modes depend on the strength of the 
vertical field, which is given by fc cr j t = 1.73Q./Va z (Balbus & Hawley 1991). The maximum growth rate of 
nonaxisymmetric modes (m ^ 0) is 0.75O at m = m max = 0.97H/uja- The cutoff wavenumber is determined 
by the strength of the toroidal field and given by m cr i t = 1.73SI/u)a (see Paper I). 
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The stable stratification acts on the MRI as the negative buoyancy and reduces its growth rate. The 
amplitude of the restoring force is denoted by the size of Nx e - The stabilizing effect due to the stratifi- 
cation on the axisymmetric (m = 0) and nonaxisymmetric (m = m max ) MRI are depicted in Figure 2. 
The normalized growth rate is given by a function of the vertical wavenumber for the cases with different 
stabilizing parameters Nxe/ft = 0,0.5, 1 and 2. The horizontal axes are measured by the logarithmic scale. 
As seen from these figures, the growth rate of MRI decreases as Nxe/Q increases. When the stabilizing 
parameter exceeds a critical value (Nxe/^ ^ 2), the growth of MRI is completely suppressed. These results 
are generally approved in spite of the radial wavenumbers. 



3.2. Effects of Quadruple-Diffusions 

Our main interest in this paper is to understand the effects of multiple diffusions on the growth of MRI. 
In general, the viscous, heat, and chemical diffusions have a major impact on the growth of hydorodynamic 
and magnetohydrodynamic instabilities (Bruenn & Dinneva 1996; Miralles et al. 2000, 2002; Bruenn et 
al. 2004). The ohmic dissipation could also affect the growth of MRI (Jin 1996; Sano & Miyama 1999). 
We examine the single and multi-diffusive effects on the growth of MRI. We restrict our discussion to the 
case with zero radial wavenumber (k r — 0). For further simplicity, we fix the azimuthal wavenumber as 
m = m max , which corresponds to the fastest growing mode of nonaxisymmetric MRI (see § 3.1). 

The dispersion relation is characterized by the normalized diffusivities defined as 

where Pe is the Peclet number, Re is the Reynolds number, and Reu is the magnetic Reynolds number. 
Here, we select the vertical Alfven speed as the typical velocity and Va z /Q as the typical lengthscale. The 
effect of the chemical diffusion is quite similar to that of the heat diffusion, and thus it is ignored for a while. 
We now consider the case with a strong toroidal field ua/ojaz = 10 4 . The characteristic length scale of 
nonaxisymmetric MRI is thus much longer than that of axisymmetric MRI. Therefore, axisymmetric MRI 
is affected by the diffusion processes more severely than nonaxisymmetric MRI. 



3.2.1. Single Diffusive Effects 

It would be useful to describe the effects of each diffusions on the MRI separately. Figures 3a, 3b and 
3c show the dispersion relation taking account of the heat, viscous, and magnetic diffusivities respectively. 
The horizontal axes of these figures are vertical wavenumber in the logarithmic scale. 

Figure 3a demonstrates the effect of the heat diffusion. The other diffusive coefficients are assumed 
to be v = r] = 0. The heat diffusion can reduce the stabilizing effect due to the stratification and revive 
the unstable growth of MRI. As described in § 3.1, the growth of MRI is suppressed by the restoring force 
of negative buoyancy. However, even with a strong stratification (TVre/^ ^ 2), MRI can grow if the heat 
diffusion is efficient. This is because the heat diffusion induces the heat exchange between the fluid element 
and the surroundings, which eliminates the entropy difference between them. Thus it reduces the stabilizing 
effects caused by the entropy gradients and the growth of MRI is revived. 

The typical wavelength X x at which the heat diffusion can affect the buoyant oscillation is estimated 
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from a simple relation, \k 2 ~ Nx e , that is 



X x = 2, x l^-. (35) 



The heat diffusion affects the modes with shorter wavelengths than A x . On the other hand, the cutoff wave- 
length of MRI is given by Amri ~ 2~kVaz/Q- Magnetic tensions stabilize the modes with longer wavelengths 
than Amri- Thus we can derive the range of poloidal wavelength which involve fast growing MRI modes in 
the case with heat diffusion and stable stratification, that is Amri ^ A < A x , in other words, 



^^W^k' (36) 

where A = A/ Amri- Here we define the fast growing MRI modes as the modes whose growth rate is of the 
order of qQ. For the case with Nt c /Q = 2, the condition for reviving the MRI is P e < 0.5. This is consistent 
with the result shown in Figure 3a. In the strongly stratified regions (iV^/fl ^> 1), such as the outer layer 
of PNSs and the radiative core of the Sun, the condition P e -C 1 is needed for the growth of MRI. 

Figure 3b shows the viscous damping of the growth of MRI. Here Nxe is assumed to be zero and the 
heat and magnetic diffusivities are neglected (x = f] = 0). The effect of the kinetic viscosity is characterized 
by the Reynolds number Re. The Reynolds number is unity when the viscous damping rate is comparable 
to the growth rate of MRI, vk 2 ~ ft. The growth of MRI is suppressed by the viscosity when Re < 1 as 
shown in Figure 3b. The critical wavelength is approximately given by 

A ^Am^ = \/^ = ^' (37) 

where A„ = v/Vt. We notice that the growth of MRI is suppressed but never stabilized by the viscous 
damping. This means that the stability condition of MRI is not affected by the viscous dissipation, which 
is discussed in § 2.3. Thus we should consider both the growth rate and the stability condition when we 
discuss the characteristics of MRI in multi-diffusive systems. 

Ohmic dissipative effects are shown in Figure 3c. The conditions N Te = and \ = v = are assumed 
in this figure. As seen from this figure, the effect of ohmic dissipation is similar to that of viscous dissipation. 
However, the ohmic dissipation can stabilize the MRI when the following condition is satisfied, 



~ v Amri V R^m ' ^ ^ 



where A^ = 2iiyJri/VL. Thus the stability condition of MRI depends on the magnetic diffusivity (see,§ 2.3). 
The ohmic dissipation becomes effective when Reu < 1 is satisfied. This is also consistent with the Figure 3c. 



3.2.2. Double Diffusive Effects 

Multi-diffusive processes often affect the stability of astrophysical objects. We investigate the double 
diffusive effect on the growth of MRI focusing on typical two cases. One is the double diffusive system with 
X 3> v (PNS-type system), and the other is that with \ 3> i] (SOLAR-type system). As described later, 
heat and viscous diffusions caused by the neutrino radiation dominate over the magnetic diffusion in PNSs 
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(x ^ v rj). On the other hand, the condition \ 3> V S> v is satisfied in the solar radiative zone. We 
restrict our discussion to the cases with zero radial wavenumber and the azimuthal wavenumber m — m max . 

Figure 4a shows the double diffusive effect on the growth of MRI in the PNS-type system. Vertical and 
horizontal axes are the same as Figure 3. The thick dashed-line shows the growth rate in the diffusionless and 
isentropic limits. The stabilization of MRI due to the stable stratification is depicted by the thick dotteddine 
(iV^/Sl = 2). The solid line shows the double diffusive effect of heat and viscous diffusions in the PNS-type 
system. We consider the case with \ ^ v i so that Pe = 10 -4 and Re = 0.01 are assumed here. Then the 
critical wavelengths are A x = 10 2 and X u = 10, respectively. 

From this figure, it is found that the stabilizing effect is reduced and the growth of MRI is revived 
for the modes with the poloidal wavelengths shorter than X x . On the other hand, the unstable growth is 
suppressed by the viscous damping for the modes with the poloidal wavelengths shorter than A„ . Therefore, 
in the PNS-type system, fast growing modes of MRI arc limited to the following range; 

max(A„, A M ri) < A < A x . (39) 

We show the relation between single and double diffusive cases. The thin dotted-line shows the single 
diffuisive case with the heat diffusion of Pe = 10~ 4 and Nt c = 2.0. The case with only the viscous diffusion 
of Re = 0.01 and Nx e — is described by the thin dashed-line. We find that the double diffusive case is 
obviously superposed by each single diffusive cases. 

Figure 4b shows the double diffusive effect on the growth of MRI in the SOLAR-type system. Thick 
dashed and dotted-lines are the same as in Figure 4a. The solid line shows the double diffusive effect in 
the SOLAR-type system (\ S> r/). When we assume Pe = 10~ 4 and Reu — 0.01, corresponding critical 
wavelengths are A x = 10 2 and A,, = 10, respectively. As seen from this figure, the unstable growth of the 
MRI are revived by the effect of heat diffusion the same as in the PNS-type system. However, the ohmic 
dissipation stablize the unstable growth of MRI in the regime A < A,,. Therefore, in the SOLAR-type system, 
fast growing modes of MRI arc limited to the following range, 

max (A,,, A M ri) < A < A x . (40) 

The SOLAR-type double diffusive case is also superposed by each single diffusive cases. The thin dashed-line 
shows the single diffuisive case with the thermal diffusion of Pe = 10~ 4 and Nt c = 2.0. The case with only 
the ohmic dissipation of Reu = 0.01 and Nx e = is depicted by the thin dotted-line. It is found that the 
condition of faster growing modes are determined by single diffusive behaviors of the disprsion relation. 

Finally, we consider the general multi-diffusive systems (x ^ 0, £ ^ 0, v ^ 0, r\ 7^ 0). It is indicated from 
the double diffusive systems that the lower limit of the wavelength for fast growing modes is determined 
by the two dissipation processes, that is the viscous and ohmic dissipation. In addition, its upper limit is 
determined by the two relaxation processes, that is the heat and lepton diffusions. Thus, in the multi-diffusive 
system, fast growing modes of MRI are limited to the following range: 

max(A„, \ rn A M R.i) < A < min (A x , A 5 ) , (41) 

where Aj = 27r x /^/iVi e . In stably stratified regions (NtcNlc > Q) such as the envelope of PNSs and 
the solar radiative core, most of unstable MRI modes are drastically suppressed by the viscous damping, 
ohmic dissipation or negative buoyancy. Only the modes satisfying the condition given by equation (41) can 
grow significantly. The properties of the systems are determined by the fastest growing modes. It is, thus, 
important to know the typical lengthscale defined by the condition (41). 
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3.2.3. Difference in Axisymmetric and Nonaxisymmetric MRI 

We have investigated the diffusive effects on the nonaxisymmetric MRI. Generally, fast growing modes 
of axisymmetric MRI are also limited to the range defined by equation (41). However, the growth rate can 
differ among the axisymmetric and nonaxisymmetric MRI in multi-diffusive cases. 

As seen from Figure 2, the growth rate of axisymmetric MRI decrease with increasing the poloidal 
wavelength, but that of nonaxisymmtric MRI remains to be constant even at k z < Q,/Vaz- Thus, in multi- 
diffusive systems, the growth of nonaxisymmetric MRI dominates over that of axisymmetric one when the 
condition, max(A„, A^) 3> Amri is satisfied. In this case, fast growing modes arc limited to the "window A" 
in Figures 2a and 2b. In the case satisfying the condition, max(A;/,A^) < Amri, on the other hand, fast 
growing modes of axisymmetric and nonaxisymmetric MRI are limited to " window £>" in Figures 2a and 
2b, and their growth rates are comparable. Therefore, it is found that characters of fast growing modes of 
axisymmetric and nonaxisymmetric MRI are determined by the strength of vertical magnetic fields when 
the sizes of diffusive coefficients are fixed. 

Note that these properties appear under the situation that typical wavelength of nonaxisymmetric MRI 
is much longer than that of axisymmetric one. The toroidal magnetic fields are generally considerd to be 
dominant over the poloidal one in the stellar context (Heger et al.2005). Therefore, these properties would 
be important in the envelope of PNSs and solar radiative core. 

4. APPLICATION TO PROTO-NEUTRON STARS 

In paper I, we have investigated the stability for the nonaxisymmetric MRI ignoring any diffusion 
processes and revealed that almost all regions of PNSs are unstable. However, we must consider the diffusive 
effects caused by the neutrino radiation in PNSs. The ohmic dissipation of hot nuclear matters can be also 
efficient in PNSs. Therefore, in this paper, the dispersion equation including multiple diffusion is derived. 
We apply it to the interiors of PNSs and reexamine the stability for the MRI. In the following, we restrict 
our discussion to the stably stratified regions of PNSs, which locates below the neutrinosphere (Janka & 
Muller 1996; Thompson & Murray 2000; Buras et al. 2003). 

4.1. Physical Quantities in Proto-Neutron Stars 

4- 1.1. Dynamical Quantities 

During the core collapse, a significant amount of differential rotation can be generated in PNSs, even if 
the precollapse core rotate rigidly (Heger et al. 2000, 2005). Numerical studies of the rotating core-collapse 
suggest that the envelope of PNSs would rotate at the angular velocity ~ 100-1000sec -1 with shear 
rate q < 1 (Buras et al.2003; Kotake et al.2003; Villain et al.2004). These values strongly depend on the 
precollapse core models. Observations of young isolated pulsars associated with supernova remnants indicate 
that they rotate with <~ 100sec -1 . Thus, we addopt fl = 100sec _1 and q = 1 as typical rotational parameters. 

From the stellar evolution calculations, we can estimate the strength of magnetic fields in PNSs. Heger et 
al. (2005) study the evolution of a magnetized massive star with rotation, and show that the precollapse iron 
core has spatially homogenious megnetic fields. Toroidal magnetic components are amplified from poloidal 
one via various dynamo processes in the stellar evolutionary phases (Spruit 1999, 2002). The toroidal 
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component is about ~ 10 9 G, which dominates over the poloidal one ~ 10 6 G. Assuming the conservation of 
magnetic flux (B oc p 2 / 3 ), the envelope of the PNSs have B = 10 13 G and B z = 10 9 G. This would be the 
lower limit of the field strengths, because the magnetic fields could be amplified by turbulent dynamos or field 
wrapping during the core-collapse (Sawai et al.2005). Typical radius and density of PNSs are R ~ 3 x 10 6 cm 
and p ~ 10 12 g cm~ 3 . The corresponding Alfvcn frequencies loa and u)a z are given by 

"^ = L0 (kItgX sxioO ' do^gV-O 3ec "'' <42) 

= ^i^i^^y'iwwy 1 ' 2 8 ° c " ■ <43 » 

The buoyancy frequencies in the stably stratified layer of PNSs are sensitive to the microscopic physics 
such as the equation of state and leptonic fraction (Buras et al.2003; Thompson et al.2005; Dessart et 
al.2005). In general, the sizes of the thermal and leptonic buoyancy frequencies are much larger than the 
angular velocity (N^ <~ N 2 » £1 2 ) and the both are of the same order below neutrinosphere (Janka & 
Miiller 1996). Thus, we assume that to be Nt ~ Nl — 10O in this section. 



If.. 1.2. Quadruple Diffusivities 

We estimate the various diffusion coefficients in PNSs. Under the cquiriblium diffusion approximation, 
we can derive the heat and chemical diffusivities analytically as shown in § 2. As the opacity Ko, we adopt 
the simplified form derived by Janka (2001); 



{m e c 2 f 

P \( k B T \2 



m u (m e c z y 



= '^Ufs^Xrasv) [cm " 1 ' <44) 

Here <r = 1.76 x 10~ 44 cm 2 is the characteristic weak interaction cross section, m e c 2 = 0.511 MeV is the 
rest-mass energy of the electron, m u is the atomic mass unit, f(a c ,Y n ,Y p ) w 1 is a function of the vector 
coupling constant a c , and the number fractions of free neutrons, Y n , and protons, Y p . The contributions 
from the neutral-current scattering and the chargcd-currcnt absorption of the electron neutrino and the anti- 
electron neutrino are included as the opacity source, but neutrino annihilation is not since the cross section 
is vanishingly small inside PNSs. Considering such the opacity, we can estimate the heat and chemical 
diffusivities as follows, 

* - 10, 1lo^)" 2 (low) 2 (i^v)"* t^-'). « 45 » 

The viscous stress is induced by the neutrinos in a hot nuclear matter, and it acts on the larger scale 
compared to the neutrino mean free path. When we assume the nondegenerate electron neutrino and anti- 
electron neutrino in degenerate nuclear matter, the heat diffusivity \ is related to the neutrino viscosity 
as 

v ~ C -fx = m 2 r / 3 p-y^[f(Y p )}-\ 

10 io/ P \- 8/3 / \ 2 f k B T ^ 



. 10 12 g cm -3 



) (lOMevJ (iMev) [cm SeC ] (4?) 
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(van den Horn & van Weert 1984; Thompson & Duncan 1993), where c p w p 1 / 3 rrh/ 3 T[f(Y p )]~ 1 is the 
specific heat at constant pressure per volume, f(Y p ) ~ 1 is a specific function of the protons number fraction 
(Thompson & Duncan 1993). This is because both of the momentum and heat are transported by neutrinos 
dominantly. 

The electric charge is transported via degenerate relativistic electrons there. Momentum transfer rate 
is mainly determined by electron-proton scattering in the case not exceeding the standard nuclear density 
~ 10 14 g cm -3 . When we consider the electron-proton collision, the electrical resistivity of hot nuclear matter 
is relatively low, K = 6x 10 -45 T 2 sec (Yakovlev & Shalybkov 1991). Therefore, the magnetic diffusivity of 
the hot nucler matters is given by 



We discuss the characteristics of the axisymmetric and nonaxisymmetric MRI in the stably stratified 
regions of PNSs. In what follows, we use the physical quantities derived in § 4.1; q = 1, O = 100 sec -1 , 
uja = 1-0 sec -1 , loaz = 10 -4 sec -1 , and = = 10 3 sec -1 . As the diffusivities, we addopt the values 
of x = 10 12 cm 2 scc -1 , £ = 3 x 10 11 cm 2 sec -1 , v = 10 10 cm 2 sec -1 , and ij = 10 -3 cm 2 sec -1 . The features of 
MRI depend on the stellar merdian altitudes. Thus, we define the polar angle 9 = tan -1 (r/2;) which denote 
the direction of the entropy and leptonic gradients. Using a polar angle 9, the radial and vertical components 
of the buoyancy frequencies are written as N r = Nsin9 and N z = N cos 9. 



First of all, we discuss the growth of MRI at the equatorial plane of PNSs (9 = tt/2) in order to show the 
diffusive effects clearly As described in Paper I, the equatorial region is stable to the MRI in the diffusionless 
case. This is because the stabilizing effect of the negative buoyancy becomes the maximum at the equatorial 
plane (Nxe, Nlc S> ft). Considering the diffusion processes, the stabilizing effect is reduced and the growth 
of MRI could be revived even at the equatorial regions. 

Figure 5a shows the maximum growth rate of MRI (k r = 0) at the equatorial plane of PNSs as a fuction 
of the azimuthal and vertical wavenumbers. The vertical wavenumber is normalized by the size of the system 
R (= 3 x 10 6 cm). In the Figure 5b, the slices of Figure 5a at m — and m = m max are depicted. From 
these figures, we find that the maximum growth rate of nonaxisymmetric MRI (~ 10 sec -1 ) is significantly 
larger than that of axisymmetric one (~ 0.01 sec -1 ). 

The range of poloidal wavelength which involve the fast growing modes of MRI is determined by the 
condition given by equation (41). Substituting the physical quantities, we obtatin the upper and lower 
bounds of the wavelength for the fast growing modes in PNSs as follows; 




(48) 



The magnetic diffusivity is extremely smaller than the other diffusion coefficients in PNSs. 



4.2. Stability of Proto-Neutron Stars 



4-2.1. Stability of the equatorial region 



10 



l 4 <A<10 5 [cm]. 



(49) 



This is surely consistent with the Figure 5. The lower bound of this range is determined by the wavelength 
for the neutrino viscosity being effective. This means that the critical wavelength of axisymmetric MRI Amri 
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is smaller than \ v . In this case, fast growing modes of MRI are limited to the "window A" in Figure 2. 
Therefore the growth of nonaxisymmetric MRI dominates over that of axisymmetric one (see § 3.2.3). 

4-2.2. Dependence of the maximum growth rate of MRI on stellar merdian altitudes 

The dependence of the maximum growth rate of MRI on stellar meridian altitudes 9 is described in 
Figure 6. Upper and lower panels show the growth rate of axisymmetric and nonaxisymmetric MRI for 
the case with different field strengths B z = 10 9 , 10 10 , and 10 12 G. As seen from these figures, the maximum 
growth rate decrease with increasing the polar angle 9. This is because the stabilizing effect due to the stable 
stratification increases with polar angle 9 (see Paper I). In addition, we can find that the maximum growth 
rate of axisymmetric MRI is strongly dependent on the strength of poloidal magnetic fields. When the 
strength of poloidal magnetic fields is weak (~ 10 9_10 G), nonaxisymmetric modes of MRI grow dominantly 
at any 9. As the strength of poloidal fields increases, the maximum growth rate of axisymmetric MRI 
increases, and then becomes comparable to that of nonaxisymmetric MRI at ~ 10 12 G. 

For the case with weak poloidal fields (B z <~ 10 9 ~ 10 G), the typical wavelength of axisymmetric MRI 
Amri is much smaller than \ v . Thus the most growing mode of axisymmetric MRI is not included in the 
range defined by equation (49) . This corresponds to the case that the fast growing modes are limited to the 
"window A" in Figure 2. Therefore, the nonaxisymmetric modes of MRI grow dominantly in this case. As 
the strength of poloidal fields increases, the length of Amri approachs to A„. This corresponds to the picture 
that the possible range of the fast growing modes shifts from "window A" to "window B" in Figure 2. Thus 
the maximum growth rate of axisymmetric MRI increases when the strength of poloidal fields increases. 

For the fast growth of axisymmetric MRI as much as that of nonaxisymmetric one, it is necessary to 
satisfy the condition \ v < Amri- Assuming the critical strength of poloidal fields satisfying this condition as 
B c , we can derive the lower limit of poloidal fields for the fast growth of axisymmetric MRI; 

B > B ^ W v- 1 „,( I5I ^)-(- 5 ^_ r )-(^) [G1 . 

Here we use the specific form of neutrino viscosity defined by equation (47). The growth of axisymmetric 
MRI is drastically suppressed by the neutrino viscosity when the strength of the poloidal magnetic fields is 
weaker than ~ 10 12 G. This critical value corresponds to that derived in Figure 5a numerically. 

Note that, in the situation with realistic quadruple diffusions, the modes with zero radial wavenumbcr 
(k r = 0) are the fastest growing modes. This is a remarkable difference from the diffusionless case that 
studied in Paper I, in which the modes with k r /k z — N r /N z are most unstable one. That is, unstable fluid 
elements can move against the stratification in neutrino loaded PNSs. This character suggests that the MRI 
itself could operate as the process leading to overturnig fluid motions (Kato 1992; Balbus & Hawley 1994). 

As a result, we conclude that stably stratified envelope of PNSs are always unstable to the MRI. Espe- 
cially, the maximum growth rate of nonaxisymmetric MRI is significantly larger than that of axisymmetric 
MRI unless the poloidal magnetic components are extremely strong (> 10 12 G). Typical growth time of MRI 
is about 100 msec. This is much shorter than the neutrino cooling time of PNSs (~ 10 sec), and thus, it 
is expected that the nonlinear growth of MRI affect the neutrino luminosity via effective heating processes 
caused by overturnig fluid motions and/or the other magnetic processes. 
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5. DISCUSSION 

The nonlinear evolution of the MRI could amplify the magnetic fields and drive MHD turbulence in 
PNSs. On the other hand, the growth of other instabilities is also expected (Bruenn & Dinneva 1996; Miralles 
ct al. 2000, 2002; Bruenn et al. 2004). In particular, the radiative-driven magnetoacoustic instability, which 
is refered as " neutrino bubble" , is focused as the seed process of the turbulent mixing and the mechanism 
for generating pulsar kicks (Blaes & Socrates 2003; Socrates et al. 2005). We compare the growth rates 
of the MRI and the neutrino bubble instability, and derive the condition for the MRI dominating over the 
neutrino bubble instability. 

The growth rate of neutrino bubble instability is given by 

7max ~ 3 x 10 2 (^j^) (iQis^ see- 2 ) (iMcv) (lO^cm" 3 ) S ° C ' ' (51) 

where B p is poloidal magnetic fields (Socrates et al. 2005). Considering the typical growth rate of MRI 
(~ qfl), we can derive the condition for the MRI dominating the neutrino bubble instability at the certain 
situation of T — 4 MeV, p <~ 10 12 g cm -3 , and g <~ 10 13 cm sec -2 : 

n > 3 x 10~ V 1 ^ . (52) 

This condition suggests that the MRI dominates over the neutrino bubble instability in the weakly magnetized 
systems with a rapid differential rotaion. 

From the condition (52), we can speculate the evolutionary scenario of PNSs. As PNSs should rotate 
differentially and be magnetized weakly right after core bounce, the growth of MRI exceeds that of the 
neutrino bubble instability. In the early evolutionaly phase, the nonlinear evolution of MRI initiates and 
sustains MHD turbulence. Thus, turbulent motions induced by MRI would enhance the neutrino luminosity 
sufficiently to yield explosions. MHD turbulence leads to angular momentum transport. Thus, the rotational 
configurations within PNSs evolve toward a rigid rotation. If the PNSs rotate rigidly, the neutrinosphcrc 
becomes oblate, and the neutrino luminosity could be enhanced in the polar region (Kotake et al. 2003, 2004). 
Thus, the polar region would be heated by neutrinos preferentially. 

After magnetic fields are amplified by the MRI and the rotation profile of PNSs becomes the rigid 
rotation, the growth rate of the neutrino bubble instability can be comparable to that of the MRI. At this 
evolutionary phase, the neutrino bubble instability can grow nonlinearly and drive turbulent mixing as is 
suggested in Socrates et al.(2005). Local luminosity enhancements, which preferentially occur in the regions 
of strong magnetic field, lead to a net global asymmetry in the neutrino emission and the young neutron 
star is propelled in the direction opposite to these regions. In addtion, the magnetic buoyancy instability 
may also grow in PNSs (Masada et al. 2006 in prep). To reveal the origin of neutron star kicks completely, 
we should study the nonlinear coupling effects of these instabilities in detail. 



6. SUMMARY 

In this paper, we investigate the stabiliy of differentially rotating, magnetized PNSs, including the effects 
of multiple-diffusions. We derive the local dispersion equation for both axisymmctric and nonaxisymmetric 
perturbations and apply it to the interior of PNSs. Our findings are summarized as follows. 

1. The heat and leptonic diffusions reduce the restoring force due to stable stratification and increase 
the maximum growth rate. The viscous damping suppresses the growth of short wavelength perturbations 
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and decrease the maximum growth rate. The ohmic dissipation stabilizes short wavelength perturbations. In 
the multi-diffusive systems, the fast growing modes of MRI can be seen only when the poloidal wavelength 
is 

max (A„, A„ A M ri) < A < min (A x , A € ) . 

The upper bound is determined by the lcngthscale that the restoring force due to the stable stratification 
is reduced by the heat or lcptonic diffusion. The lower bound is determined by the lengthscale for ohmic 
dissipation or viscous damping being effective. 

2. The necessary conditions of the stability for MRI in multi-diffusive case are derived from our simplified 
dispersion equation. The conditions are similar to that of Menou et al.(2004), in which only the stability for 
axisymmetric MRI is investigated. 

3. MRI can grow in realistic neutrino loaded PNSs. For the case with » B z , the growth of 
nonaxisymmctric MRI dominates over axisymmetric one unless the vertical field is extremely strong [B z > 
10 12 G). Typical growth time of fast growing MRI modes is about 100 msec, which is much shorter than the 
neutrino cooling time of PNSs. 

4. The condition that the growth of the MRI dominates over the neutrino bubble instability is given by 

fi > 3 x 10- V'Sp , 

at T <~ 4MeV, p <~ 10 12 g cm -3 and g ~ 10 13 cm sec" 2 . Right after the core bounce, the growth of the MRI 
is dominant over the neutrino bubble instability. After the strong magnetic field is generated by the MRI 
and the rotational configuration shift to a rigid shape, the growth rate of the neutrino bubble instability 
could become larger than that of the MRI. 

We thank V. Urpin and S. Nagataki for helpful discussions. TS is supported by the Grant-in-Aid 
(16740111, 17039005) from the Ministry of Education, Culture, Sports, Science, and Technology of Japan. 
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Fig. 1.— Unstable growth rate ^s(u>) in the diffusionlcss case is plotted as a function the azimuthal wavenum- 
ber m and the vertical wavenumber k z . For simplicity, we select the case of zero radial wavenumber k r = 0. 
The growth rate is normalized by the angular velocity f2, and the wavenumbers are normalized by loa/Q and 
Vaz/Q, respectively. In this figure, we assume the shear parameter q — 1.5, the effective thermal and leptonic 
buoyancy frequencies Nx e — Nl c = 0. Quadruple diffusivities are assumed to be zero (x = £ = ^ = 77 = 0). 
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Fig. 2. — Dependence of the growth rate on the stabilizing parameter iV<r e /f2 for the case of (a) axisymmetric 
MRI (m = 0), and (b) nonaxisymmetric MRI (m = m max ). The growth rate is shown as a function of the 
normalized vertical wavenumber k z VAz/Q for the case Nre/Q = 0, 0.5, 1, and 2, respectively. Note that the 
horizontal axis is measured by the logarithmic scale. Model parameters are the same as those in Figure 1. 
The growth rate of axisymmetric and nonaxisymmetric MRI decrease as the stabilizing parameter increases. 
When Ntc ^ 2f2, the growth of the MRI is completely stabilized by the stable stratification. 
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Fig. 3. — Single diffusive effects of (a) heat, (b) viscous, and (c) magnetic diffusions on the growth of MM. 
Normalized growth rates are shown as a function of the normalized vertical wavenumber k z VAz/Q- We use 
the Peclet number Pe, the Reynolds number Re and the magnetic Reynolds number Reu as the measures 
of the effects of heat, viscous, and mangetic diffusivities respectively. For simplicity, we select the branch 
with m = m max and k r = 0. This corresponds to the fastest growing branch of the nonaxisymmetric MRI. 
The other parameters are the same as in Figure 1. 
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Fig. 4. — Double diffusive effects in (a) the PNS-type systems and (b) the SOLAR-type systems. In the 
PNS-type systems, the heat and viscous diffusivities are sufficiently larger than the magnetic diffusivity 
(x ^ v ^ v)- C> n the other hand, in the SOLAR-type systems, the heat and magnetic diffusivities are 
dominant over the viscous diffusivity (x 3> 77 3> ^). Unstable growth rates are shown as a function of the 
normalized vertical wavenumber k z VAz/Q- Model parameters are the same as those in Figure 3. 
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Fig. 5. — (a) Three-dimensional plot of the unstable growth rate as a function of the azimuthal wavenumber 
m and the normalized vertical wavenumber k z R, and (b) slices of the three-dimensional plot at m = 
and m = m max . The size of the system is assumed to be R = 3 x 10 6 cm. The modes with m = are 
corresponding to the axisymmetric MRI (solid line), and the other is the nonaxisymmetric MM (dashed- 
line). We choose the model parameters as q = 1.0, SI = 100 sec -1 , Nt = Nl = 10f2 sec -1 , u>a = 1-0 sec -1 , 
and u>Az = 10 -4 sec -1 . Quadruple diffusivities are assumed as \ = 10 12 cm 2 sec -1 , ^ = 3 x 10 11 cm 2 sec -1 , 
v = 10 10 cm 2 sec -1 , and r\ — 10 -3 cm 2 sec -1 , respectively. These figures represent the characteristics of 
MRI at the equatorial plane of the PNSs. 
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Fig. 6. — The maximum growth rate of the MRI as a function of the polar angle 9 for the cases of B z = 
10 9 , 10 10 , and 10 12 G. The upper panel shows the behavior of axisymmetric MRI (m = 0), and the lower panel 
is nonaxisymmetric MRI [m ^ 0). We search the fastest growing modes solving our simplyfied dispersion 
equation (31). The strength of toroidal fields is assumed to be constant in all cases, B$ = 10 13 G. The other 
parameters are the same as those in Figure 5. Note that the vertical axis of upper figure is measured by the 
logarithmic scale, but that of bottom figure is measured by the normal scale. 



